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ABSTRACT 

The  results    of    two   different   methods    for    determining   the 
rigidity   modulus    of  a    soft    sediment    are   compared.       In    one   method 
the   resonant    characteristics    of   a    torsionally    oscillating    rod 
which   are   sensitive   to   the    shear    acoustic    impedance   of   the   sedi- 
ment   in   which   the   rod   is    imbedded   determine   the    complex   rigidity. 
The   second  method  utilizes    the   observation   of   the  phase   velocity 
of  an    interface   wave  at    the  water-sediment   boundry.      Shear    wave 
speeds    computed   from   the    experimental    data   from   both  methods    are 
quite   similar    in   magnitude.       For    the    sediment    used   here,    the 
average   value   of   shear   wave    speed  determined  from   the    interface 
wave    experiment   was    33  m/sec   while   the   shear    wave   speed   determined 
from   the   measured  rigidity   was   29   m/sec.      The   difference   lies   with- 
in   experimental    uncertainty.      Trends    in    the   mass-physical   properlies 
of   the   sediments    are    investigated  by    comparing   graphically    the    de- 
pendence  of   both    the  real    and   imaginary   parts    of    the   complex   rigid- 
ity   on    density,    porosity,    sound   speed,    silt   and   clay  percentages ., 
Poisson's    ratio   and   density    times    sound   speed  squared. 
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I.   INTRODUCTION 

Growing  concern  for  the  accurate  and  comprehensive  description 
of  the  acoustic  reflection  and  attenuation  characteristics  of  the 
ocean  bottom  is  continually  stimulated  by  attempts  to  fill  the 
large  gap  in  predictibility  of  the  actual  reflection  losses.   The 
development  of  a  model  for  the  reflection  coefficient  which  more 
accurately  describes  the  physics  of  transmission  and  reflection 
should  result  in  improved  capability  for  predicting  transmission 
loss  for  sound  in  the  ocean  in  those  cases  where  bottom  interaction 
is  important,  such  as  passive  and  active  sonars  in  shallow  water 
and  active  sonars  using  the  bottom  bounce  mode. 

Historically,  the  ocean  bottom  has  most  often  been  modeled  as 
an  interface  between  two  fluids  of  differing  densities  and  sound 
speeds.   Usually  referred  to  as  the  Rayleigh  model,  it  is  generally 
accepted  because  normal  surficial  sediments  have  high  porosities 
and  low  rigidities  which  approximate  those  of  a  fluid.   The 
Rayleigh  model  is  sometimes  found  to  be  deficient,  especially  at 
lower  grazing  angles  [3,9]. 

Some  have  accounted  for  this  disparity  by  presuming  that  the 
sediment  posseses  a  finite  rigidity  and  hence  is  capable  of  sup- 
porting a  shear  wave  [ 3] •   Using  a  viscoelastic  model  and  complex 
Lame1  constants,  one  is  capable  of  describing  the  formation  of 
shear -inter face  waves  and  absorption  of  compressive  sound  waves 
in  the  sediment.   Bucker ,  et  al.  [3]  finds  better  agreement  with 


experimental   measurements    for    the   theoretically   calculated  bottom 
loss   using   this    viscoelastic   model   and  associated  sediment   para- 
meters   obtained  by  Hamilton   [Q,10,ll]    than   for    the   Rayleigh   model 
in   a   number    of   cases. 

In   order    to  obtain   both   complex  Lame'    constants    for    the   vis- 
coelastic  model    it    is    required   that,    in   addition    to   the   measure- 
ment   of    sound  absorption,    both    shear    wave   speed  and   compressional 
wave   speed  be  measured.      Hutchins    Ll2]    demonstrated   the   feasibility 
of  measuring   the   complex   shear    modulus    of   soft   sediments    (artifi- 
cial)   by   noting   the    effects    caused  by    the   rigidity    on    the   resonant 
behavior    of   a    torsionally    oscillating    rod    imbedded   in    the   sediment. 
Utilizing   a    laboratory  prepared  kaolin-water    mixture,    he    developed 
the   viscoelastometer    technique    in   the   frequency    range   of   38.3    to 
38.9  KHz. 

Cohen   L5]»    by    improving   the   design    of    the    instrument,    increased 
the   range   of    operation    to    5.8   to   38  KHz.      Cohen   used  both    kaolin- 
water    and  bentonite-water    combinations   as   artificial   sediments. 

Utilizing   Hutchins'    and  Cohen's    observations    that    soft, 
clayey,    simulated   sediments    did,    in   fact,    have  a   measureable   shear 
modulus,    Bieda    [2]    continued  this   research  by    developing   a    new 
viscoelastometer    suitable   for    use   in   making    in    situ  measurements 
and   the   subsequent    laboratory   analysis    of   20   sediment   cores    ob- 
tained  in   the  Monterey    Bay   area. 

The   results    obtained  from  use    of   these  probes    [2,5,12]    are 
quite   encouraging  and  are   certainly   within    reasonable    bounds    for 


this  type  of  sedimentary  material  (E.L.  Hamilton,  personal 
communication).   However,  there  has  been  no  independent  deter- 
mination of  the  rigidity  for  comparison  and  verification  of  the 
method.   Consequently  it  was  required  that  experiments  be  con- 
ducted that  would  test  the  validity  of  the  probe  method.   Suc- 
cessful development  of  such  a  probe  would  make  the  necessary 
measurements  of  ocean  bottom  parameters  a  more  convenient 
operation . 

Another  method  for  determining  the  elastic  constants  involves 
measuring  directly  both  the  compress ional  wave  and  shear  wave 
speeds;  knowing  these,  it  is  a  simple  matter  to  determine  the  real 
elastic  constants.   The  imaginary  parts  of  the  elastic  moduli  must 
be  obtained  from  measurements  of  the  attenuation  of  the  compress- 
ional  and  shear  waves.   The  former  method  is  used  in  this  experi- 
ment as  an  independent  means  for  validating  the  measurement  of 
the  real  components  of  the  rigidity  made  by  the  viscoelastometer 
developed  by  Bieda  [2].   The  method  described  is  essentially  the 
same  as  that  used  by  Hamilton,  et  al.  [ll]  in  measurement  of  the 
interface  wave  speed  and  compressional  wave  attenuation  in  sedi- 
ments in  the  San  Diego  area  . 


Some  controversy  appears  to  exist  as  to  just  what  title  should 
be  applied  to  the  interface  wave  that  exists  at  the  suface  between 
an  elastic,  solid  bottom  and  the  overlying  water.   Although  gener- 
ally referred  to  as  Stoneley  waves,  it  appears  that  Stoneley  C 17] 
never  included  the  possibility  that  the  upper  medium  could  be  a 
fluid.   Other  people  tend  to  choose  the  name  Scholte  waves  since 
Scholte  [l6 J  generalized  the  nomenclature  to  include  fluid-solid 
interface  waves.   Since  the  environment  in  which  the  author's  data 
were  taken  did  not  entirely  approach  an  elastic-solid  bottom  of 
infinite  extent  below  a  infinite  half-space  of  homogeneous  fluid, 
it  was  decided  to  use  the  term  "interface  wave"  when  describing  the 
elastic  wave  existing  at  the  interface  between  the  bottom  and  the 
overlying  water. 
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II.   EQUIPMENT  AND  PROCEDURES 

A.  LOCATION  OF  THE  EXPERIMENT 

Cohen  [5]  has  shown  the  viscoelastometer  to  be  ineffective  in 
sediments  of  large  mean  grain  size,  possibly  because  the  grain 
dimension  becomes  comparable  to  the  shear  wavelength  in  the  sandy 
sediments  and  therefore  the  assumption  of  homogeneity  becomes 
invalid.   Thus  there  was  a  need  for  finding  a  sediment  which  was 
silty  clay  or  clayey  silt. 

There  are  vast  areas  of  mud  bottoms  in  the  Monterey  Submarine 
Canyon  area  of  Monterey  Bay,  California,  that  would  normally 
suffice  for  this  type  of  experiment.   However,  since  the  services 
of  a  submersible  platform  were  not  available  the  evaluation  had 
to  be  made  at  a  quite  shallow  depths  (1  to  1.5  m)  near  shore.   A 
suitable  mud  sediment  site  was  found  in  the  tidal  flats  of  the 
Kirby  Park  area  of  Elkhorn  Slough,  a  few  miles  north  of  the  city 
of  Monterey . 

B.  INTERFACE  WAVE  MEASUREMENT  TECHNIQUE 

Direct  measurement  of  the  speed  of  shear  waves  in  very  soft 
marine  sediments  is  made  virtually  impossible  at  acoustic  fre- 
quencies because  of  the  rapid  attenuation  of  the  shear  waves  and 
the  difficulty  of  coupling  a  suitable  transducer  to  the  medium. 
It  has,  however,  been  determined  by  Strict  and  Ginzbarg  [l8J  and 
by  Hamilton,  et  al.  [ll]  that  the  ratio  of  the  shear  wave  speed 
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to  the    intetface   wave   speed  remains   relatively   constant    for   a 
given   value    of   the    ratio   of   the   two   densities    of  adjoining   media 
irrespective   of   other    varying   bottom   constants. 

While   compressional   wave  speed  and  wet    density    of   the   sediment 
can    easily    be  measured   in    the   laboratory   and   corrected  to    in   situ 
values,    it   was    required   that   the    interface  wave   speed  be   measured 
and   corrected   to   shear    speed  as   previously   accomplished   by 
Hamilton,    et   al.    [ll]. 

An    interface   wave,    generated  by    the   detonation   of   a   blasting 
cap   at    the  water -sediment    interface,    is    examined  using   three   or 
four    vertically-polarized,    velocity-sensitive   geophones   mounted  at 
known    intervals    on    the    interface.      An    aluminum   rod   and   flexible 
lines    served  as    a    convenient    support   and    the   means    for    spacing    of 
the   geophones   and   the  electrically-fired  blasting   caps   which  were 
mounted  under    a   heavy   brass   plate. 

During    the   firing    of   these   blasting    caps,    the   response   of    the 
geophones    is   monitered  at   a   remote   station.      The  signal    from   the 
geophones   passes    through   an    impedance-matching    transformer,    a    set 
of   attenuator   pads    and  a    low-noise  amplifier    capable    of   as   much   as 
50   db  gain    to  a   Honeywell   Visicorder    where    the   subsequent   amplitude 
versus    time   responses    are  plotted  as   shown    in   Figure   2.      A  sketch 
of   the    operational    assembly    is    shown    in   Figure   1. 

The  "Visicor der    traces    are   studied   to  determine   the  points    of 
constant  phase    in    each    of    the   geophone    traces    (point   a    in    Figure   2) . 
This    is    generally    best    accomplished  by   attempting    to  match   waveforms 


12 


which  is  somewhat  aided  by  irregular  spacing  of  the  geophones  and 
adjustment  of  gain  so  that  amplitudes  are  comparable.   Having 
found  the  points  of  constant  phase,  and  using  convenient  time 
reference  marks,  the  interface  wave  speed  is  easily  determined. 

Having  obtained  the  interface  wave  speed,  the  graphs  provided 
by  Strict  and  Ginzbarg  [l8]  are  entered  with  the  ratio  of  the 
density  of  the  overlying  sea-water  to  the  density  of  the  sediment, 
The  graphs  are  plotted  with  the  ratio  of  interface  wave  speed  to 
shear  wave  speed  versus  the  ratio  of  shear  wave  speed  to  compress- 
ional  wave  speed,  varying  the  parameter  ratio  of  water  density  to 
sediment  density.   Since  both  the  ordinate  and  abscissa  are  ratios 
involving  the  unknown,  shear  wave  speed,  an  additional  assumption 
about  the  magnitude  of  shear  wave  speed  allows  one  to  accurately 
determine  the  necessary  value.   As  the  ratio  of  shear  wave  speed 
to  compressional  wave  speed  approaches  0.5,  the  individual  graphs 
all  begin  to  approach  a  horizontal  asymptote  which  may  conveniently 
be  read  from  the  ordinate  with  a  maximum  possible  error  of 
5  percent. 

Having  determined  the  value  of  the  ratio  of  the  interface  wave 

speed  to  shear  wave  speed,  and  knowing  the  interface  wave  speed, 

it  is  easy  to  compute  the  shear  wave  speed  (V  =  V./A,  where  A  is 

the  ratio,  V  the  shear  wave  speed  and  V.  the  interface  wave 
'   s  r         1 

speed) . 
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C.      THE  VISCOELASTOMETER 

1 .  Equipment 

The  apparatus    designed  by    Bieda    L2J    is    a   result    of  a 
series    of    improvements    on    the   original    devices    constructed  and 
evaluated  by   Cohen   [5]    and  Hutchins    [ 12] .       It    consists    of   a 
10^- inch    long,    ?g-inch    diameter    cylindrical    rod  which    is    excited 
in   torsional    oscillation   by  a    centrally  mounted  barium   titanate 
transducer.      The   rods,    constructed   from  a    constant   modulus    alloy, 
Nl-Span-C,    were  heattreated  and   highly  polished.      The    electrical 
leads    to    the    transducer    are  passed   through   a    supporting    yoke  which 
supports    a    clamp    that    holds    the    transducer    in   position.      The   clamp 
is    secured   to   the   transducer    at    its    center,    a    node   of   torsional 
motion   for    all    resonant   modes    in    the   rod.      A    sketch   of    the 
apparatus    is    shown    in    Figure    3-       Various    compounds   were   used  to 
maintain    the   water-tight    integrity    of    the    electrical    connections, 
generally   RTV    insulating   rubber    was    deposited    in   the    large   crevices 
at    the   junction    of   the   clamp   and   the    transducer    while    the    other 
boundries    were   coated   with   varnish    or    glyptol. 

A    Dranetz    100B    Impedance-Admittance   meter ,    a    suitable 
oscillator    and  a   frequency   counter    were  used   to  measure   the   elec- 
trical  conductance  and   the  resonant    frequencies. 

2 .  Theory   and  Method   of  Measurement 

"When    the   viscoelastometer    is    excited  at    its   mechanical 
resonance,    a    torsionally    oscillating   standing   wave   exists    along 
the  axis    of    the   rods    and   transducer.      The    torsional   motion    of    the 
walls    of    the   rod  generates   a    shear    wave  which   propagates    radially 
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outward  into  the  sediment.   The  radiation  impedance  of  the  sur- 
rounding sediment  has  a  modifying  effect  on  the  magnitude  of  the 
electrical  resistance  of  the  probe  and  its  frequency  of  mechanical 
resonance,  compared  to  those  values  measured  in  air  (approximately 
no-load  conditions).   The  change  in  resistance  and  resonant  fre- 
quency are  thus  direct  functions  of  the  rigidity  of  the  sediment  Ll4] 

Under  the  assumption  that  the  shear  waves  that  propagate 
outward  from  the  viscoelastometer  have  wavelengths  much  smaller  than 
the  radius  of  curvature  of  the  rods  and  that  the  amplitude  of  these 
shear  waves  is  very  rapidly  attenuated,  the  waves  propagated  in  the 
sediment  may  be  considered  to  be  plane  waves  with  the  formula: 

T   =  T      exp    [  \|rTfp/T]    (l+i)z  ]  (1) 

where  T    is    the   shear    stress   amplitude,    T      the    initial   amplitude, 
f   the   frequency,    and  77    the   shear    viscosity.      The  harmonic    time- 
varying   function   has    been   suppressed  for    sake   of  convenience. 

In   a   Newtonian   fluid,    the   coefficient    of  shear    viscosity 
is  77    =  T/S,    where   S    is    the  shear    strain  and   S   =   bS/ftt.       In   our 
viscosity   model   a    complex  rigidity    (G    )    has    been  assumed  along 
with  a   complex   viscosity    (77    ).      Thus, 

*?c   =  Tl1    -    i772   and  Gc    =   Gx    +    iG2.  (2) 

In   a   Newtonian    fluid,   77      =   0,    and  T],    =  0    in   a   perfectly   elastic 
solid.      Since  77      =  T/i  CO  S   =   -    iG  /  CO    for    simple   harmonic   motion, 
where  CO  is    the  angular    velocity, 

Gl    =  WTfe   and  G2    =  WT?1  '  ^ 
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The  specific  radiation  impedance  (Z)  of  the  sediment  is 
the  ratio  of  shear  stress  to  particle  velocity.   Following 
Mason  [l4l  and  McSkimin  [l5J»  for  a  Newtonian  fluid 


Z  =  R  +  iX  =  ^TTf  r?p  (1  +  i)  (4) 

where  R  and  X  are  the  specific  acoustic  impedance  and  reactance 
of  the  sediment,  respectively,  and  p  is  the  sediment  density.   For 
the  case  when  viscosity  is  complex,  the  following  relations  result 


2   2 
"l  *  COO  ™2  LOO 


2   2 

~   -  R  -X  _   _  2RX 

Gl      ~  G2  ~  ~     * 


(5) 


Thus,  it  becomes  apparent  that  determination  of  the  complex 
rigidity  or  viscosity  requires  only  a  knowledge  of  real  and 
imaginary  parts  of  the  specific  acoustic  impedance  for  plane 
shear  waves . 

Mason  [l4]  has  shown  that  constants  of  proportionality  can 

2 
be  determined  which  relate  the  total  acoustic  resistance  (R+)  to 

the  change  of  electrical  resistance  ( AR  ),  i.e.,  AR   =  K, R^ ,  and 
v  K      e' '      '    e    1  t 

the  total  acoustic  reactance  (X  )  to  the  change  of  resonant  fre- 
quency (Af),  i.e.,  Af  =  K_X  ,  when  the  loading  on  the  rod  is 
changed  from  air  to  the  medium  under  study.   Here,  R  =  R  /l  and 
X  =  X  /I ,  where  1  is  the  length  of  the  probe  immersed  in  the 
sediment.   Thus,  measurement  of  the  acoustic  impedance  is  reduced 


2 

Transducer  which  drives  the  rod. 
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to  measurement  of  the  change  observed  in  the  frequency  of 

resonance  and  in  the  resonant  value  of  electrical  conductance 

when  the  probe  is  operated  in  the  material  to  be  analyzed. 

Evaluation  of  the  constants  K   and  K_  is  accomplished  by 

measurement  of  AR   and  Af  in  Newtonian  calibration  fluids.   Since 

e 

77      =   0    in   a   Newtonian    fluid,    R   =   X;    knowing    the    viscosity   and 
density   of   the    test    fluid,    the  procedure   for    calculating   K      and 
K_    is    straightforward. 
3.      Core  Analysis 

In    the   same  relative    location   of   the    interface   wave   experi- 
ment   ten   ceres    were   taken    for    purposes    of   viscoelastometer    analysis, 
density   and   compressional   wave   speed  measurement.      While   maintaining 
the    core  at    its    original    temperature,    the   complex   rigidity    is    de- 
termined using    Bieda's    equipment.      The   measured   real   part    of   rigid- 
ity     and   the   observed   density   are   used   to   determine   a    shear    wave 
speed,    i.e.,    Vg    =    (G^p)2. 

The   core   samples    were  analyzed  at   two   of    the   resonant 
frequencies,    27    kHz   and    38    kHz,    due    to  anomalous    behavior    at    the 
lower    frequencies,    especially    in    the    stiffer    samples.      Since   the 
constants   K..    and  K„    tend   to   vary   with    frequency    and   viscosity,    an 
average   value   was   used  at    each   frequency   from  calibrations    in 
fluids   with   viscosities    roughly    equivalent   to    those   expected   in 
the   sediments.      Table    III    contains    examples    of   the   variance   of 
these   constants    over    temperature   and   viscosity. 


17 


The  computed  shear  velocities  are  averaged  and  are  com- 
pared to  the  value  for  shear  velocity  obtained  in  the  interface 
wave  experiment  in  Table  II. 

D.   COMPRESSIONAL  WAVE  SPEED  MEASUREMENT 

Compress ional  wave  speed  is  measured  in  each  of  the  sample 
cores  by  determination  of  the  time  required  for  a  2  MHz  com- 
press ional  pulse  to  pass  through  a  known  length  of  sediment. 
Barium  titanate  disk  transducers  are  used  as  a  transmitter  and 
as  a  receiver ,  while  the  transmitted  pulse  width  and  frequency 
are  varied  in  close  proximity  to  2  MHz  in  order  to  obtain  resonant 
operation  of  the  transducers. 

The  two  transducers  are  mounted  facing  each  other  as  shown 
in  Figure  4.   The  sample,  approximately  3  inches  long,  is  inserted 
in  the  holder  which  is  in  contact  with  one  transducer;  the  other 
transducer  can  then  be  moved  to  contact  the  other  end  of  the 
sample.   A  dial  indicator  is  used  to  measure  the  separation  while 
a  time  mark  generator  and  the  output  of  a  triggered  pulsed  oscillator 
are  displayed  on  the  face  of  a  dual  trace  oscilloscope.   The  tem- 
perature is  maintained  at  a  constant  value,  identical  to  that 
measured  during  the  interface  wave  experiment. 

Calibration  of  the  velocimeter  is  accomplished  by  measuring 
a  sequence  of  time  delays  for  various  transducer  spacings  in  a 
mixture  of  22  parts  ethanol  to  100  parts  distilled  water.   This 
combination  was  shown  to  be  relatively  insensitive  to  small  tem- 
perature changes  in  the  neighborhood  of  the  temperatures  of  this 
experiment  [2l]. 
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The  sediment  sound  speeds  and  the  speed  of  sound  in  the  water 
samples  (computed  using  Wilson's  equation)  taken  in  the  vicinity 
of  the  experiments  are  shown  in  Table  I. 


E.   WET  DENSITY  AND  POROSITY 

Wet  density  is  measured  by  weighing  a  known  volume  of  the 
saturated  sediment.   The  weight  of  this  same  sample  after  heating 
in  an  oven  at  105  C  for  a  period  of  24  hours  is  again  determined. 
If  it  is  assumed  that  the  weight  lost  is  completely  derived  from 
evaporation  of  the  entrapped  water  and  that  the  weight  of  1  mX  of 
water  is  lg,  the  porosity  (P)  is  easily  calculated, 


P  = 


Volume  of  Voids 
Total  Volume 


(6) 


F»   CALCULATED  PARAMETERS 

The  shear  wave  speed  is  calculated  by  use  of  the  real  part  of 
the  measured  rigidity,  G..  ,  and  the  sediment  density,  p,  in  the 
formula: 


v|f  ' 


(7) 


shear    wave   speed   is    also  computed   from   inteface  wave   spee"d  by    the 
formula: 


Vs    =  V./A    , 


(8) 


where  A  is  the  constant  of  proportionality  derived  from  the  graphs 
of  Strict  and  Ginzbarg  Ll8j.  Poisson's  ratio  0",  is  calculated  by 
the   relation: 


o-\ 


1  - 


«v/V.)--i) 


(9) 
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where  V      is    the   compress ional   wave   speed  and  V      is    the   shear    wave 
P  s 

speed.      The   measured  and   computed  physical    constants    are   shown 
in  Table   I . 

G.       SAND- SILT -CLAY    RATIOS:    MINERALOGY 

Sand-silt-clay   ratios    were   calculated  at   the  Naval  Postgraduate 
School    by    the    sieve  pipette   method.      Clay   mineralogy   analyses, 
performed  by    Dr.    W.R.    Bryant,   Texas   A&M   University,    are   shown    in 
the  Appendix. 
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III.   LIMITATION  OF  THE  METHODS 

i 

A.   RESONANT  VISCOELASTOMETER  TECHNIQUE 

While  usage  of  this  probe  is  not  strictly  limited  to  the 
laboratory,  a  very  strong  dependence  on  temperature  of  the  con- 
stants K  and  K„  (Table  III)  requires  that  measurements  for  com- 
parison be  made  in  a  temperature-controlled  environment.   To 
allow  for  correct  measurements  to  be  made  from  2-3/4  inch  diameter 
cores,  the  assumption  is  made  that  the  shear  wave  length  is  much 
smaller  than  the  distance  from  the  probe  to  the  core  liner.   It 
was  noted  that  at  the  lower  frequency,  especially  in  stiff er 
sediments,  erratic  readings  were  generally  obtained,  possibly  from 
the  violation  of  this  assumption.   This  limitation  should  not  pose 
any  great  problems  when  the  probe  is  used  in  situ  as  there  should 
then  be  no  surrounding  barrier  in  coupling  to  the  motion  of  the 
oscillating  rods. 

It  was  also  noted  that  the  values  of  K,  and  K„  vary  with 
changing  viscosity  as  well  as  frequency  and  temperature.   Con- 
sequently, the  final  constants  used  were  obtained  from  an  average 
of  those  values  obtained  over  the  lower  values  of  viscosity 
(1-15  poise).   This  could  cause  an  error  in  K..  of  no  more  than 
30  percent  and  perhaps  110  percent  in  the  value  of  K„ . 

Another  obvious  limitation  to  the  laboratory  analysis  is  that 
no  matter  how  carefully  a  core  is  taken,  there  is  always  some 
structural  distortion  of  the  sample  as  well  as  some  drying  before 
analysis  can  be  made. 


21 


B.   INTERFACE  WAVE  TECHNIQUE 

The  most  serious  breach  of  assumptions  occurs  in  the  measure- 
ment of  the  interface  waves  at  shallow  depths.   The  determination 
of  the  interface  wave  speeds  are  all  predicted  on  an  interface 
between  liquid  and  elastic-solid  half  spaces  L4»7j-   In  the  mea- 
surement of  the  lower  frequency  (about  10  Hz)  interface  waves  in 
shallow  water  it  would  seem  that  the  effect  of  the  overlying  water 
is  quite  small,  since  the  wavelength  for  10  Hz  compr essional  waves 
in  water  is  so  large  compared  to  the  depth.   Thus,  the  situation 
is  very  nearly  that  normally  encountered  for  Rayleigh  waves.   This 
effect  could  cause  the  ratio  of  the  interface  wave  speed  to  the 
shear  wave  speed,  A,  to  approach  a  value  of  0.92  L4 J •   This  could 
yield  a  possible  error  in  the  measured  shear  wave  speed  of  as  much 
as  10  percent. 

Higher  frequency  compressional  waves,  although  much  faster, 
occasionally  tend  to  mask  the  arrival  of  the  interface  wave  due  to 
the  exponential  tail  of  the  explosion.   Often  this  masking  caused 
some  uncertainty  as  to  the  actual  arrival  time  of  the  interface 
wave.   The  standard  deviation  of  the  measured  values  of  interface 
wave  speed  are  about  11  percent  of  the  average  measured  value. 

Since  the  distance  between  the  geophones  enters  directly  into 
the  computation  of  the  interface  wave  speed,  the  uncertainty  of 
about  2 "cm  could  yield  an  error  of  about  5  percent  in  the  result. 

Since  only  quantitative  information  describing  the  relation- 
ship between  and  interface  wave  an   a  solid  shear  wave  is  that 
given  by  Strict  and  Ginzbarg  [l8],  the  author  was  forced  to  make 
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the  approximation  that  shear  wave  speed  would  always  be  less  than 
one-half  the  compressional  wave  speed  in  order  to  obtain  a  value 
for  the  needed  ratio,  A,  as  noted  in  section  II-B.   This  uncer- 
tainty in  the  constant,  A,  could  yield  an  error  of  no  greater 
magnitude  than  5  percent. 
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IV.   RESULTS  AND  DISCUSSION 

Table  II  presents  a  comparison  of  the  two  methods  of  measuring 
rigidity  by  showing  the  average  values  of  the  shear  wave  speed  as 
calculated  from  measurements  of  the  real  part  of  the  rigidity  and 
the  shear  wave  speed  computed  from  the  interface  wave  experiment. 
Extremely  good  correlation  between  the  two  values  tends  to  confirm 
the  validity  of  the  torsional  viscoelastometer  results.   It 
must  be  pointed  out,  however,  that  the  viscoelastometer  was 
operated  at  frequencies  of  27  kHz  and  38  kHz  while  the  interface 
wave  frequencies  were  about  10  Hz.   Although  there  is  a  great 
disparity  in  frequencies,  there  is  no  reason  to  expect  that 
dispersion  of  shear  speed  should  occur  LllJ.   The  viscoelastometer 
technique  is  subject  to  a  maximum  precision  of  30  percent  in 
calculating  the  real  part  of  the  rigidity,  largely  arising  from 
an  inability  to  accurately  measure  all  the  parameters.   While 
these  uncertainties  of  the  method  point  toward  the  need  for 
stricter  maintenance  of  conditions  during  viscoelastometer  mea- 
surements and  during  calibration,  they  are  not  of  sufficient 
magnitude  to  discredit  the  method. 

The  mechanical  disturbance  in  coring  and  in  transporting  the 
samples. very  likely  disturbs  the  physical  properties  by  an  inde- 
terminant  amount.   It  is  also  known  from  a  visual  inspection  of 
the  cores  that  the  bottom  in  this  area  was  far  from  being  homo- 
geneous, in  either  the  horizontal  or  the  vertical  plane. 
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The  measured  physical  parameters  for  the  ten  cores  taken  tend 
to  follow  the  same  general  patterns  (Figures  5  to  19)  as  did  those 
presented  by  Bieda  [2]  with  the  major  exception  that  the  cores 
taken  in  the  slough  possess    less  than  7  percent  sand,  making 
them  much  more  acceptable  for  viscoelastometer  analysis. 

It  is  significant  to  note  that  the  rigidities  obtained  ranged 

between  5-1  x  10   to  2.7  x  10   dynes/cm  for  the  real  part  and 

6  6         ? 

1.1  x  10   to  7.9  x  10   dynes/cm   for  the  imaginary  part.   These 

values,  as  pointed  out  by  Bieda  [2],  are  significantly  lower  than 

those  found  by  Hamilton  [ll],  but  additional  evidence  continues 

to  point  toward  the  direction  of  lower  rigidities  in  surficial 

sediments  [l,6j. 

The  trends  noted  by  Bieda  were  also  observed  in  this  research; 

the  real  part  of  rigidity  increases  with  increasing  density, 

2     ... 
increasing  percent  silt  and  with  increasing  pv   .   Rigidity  de- 
creases with  increasing  porosity  and  percent  clay.   The  same 
trends  are  noticeable  but  less  marked  with  the  imaginary  part  due 
to  the  greater  degree  of  scatter  caused  by  experimental  imprecision 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  relatively  close  agreement  of  the  two  values  of  shear 
wave  speed  strongly  tends  to  validate  the  usefulness  of  the 
viscoelastometer  technique  for  measuring  the  real  part  of  the 
complex  rigidity  in  soft,  fine  grained  surficial  marine  sediments. 

The  mass-physical  trends  noted  by  Bieda  [2]  are  seemingly 
substantiated  by  these  observations  but  require  a  great  deal  more 
study  to  completely  classify  definite  behavioral  functions. 

In  order  to  further  exploit  the  use  of  the  viscoelastometer, 
it  is  recommended  that  further  efforts  be  applied  to  the  design 
of  the  viscoelastometer  in  order  to  minimize  the  effects  on  its 
calibration  of  grain  size,      stiffness   and  temperature. 

Further  effort  is  also  indicated  toward  improving  the  precision 
in  measuring  the  imaginary  part  of  the  complex  rigidity  and  toward 
its  validation,  utilizing  a  measurement  of  the  sound  absorption  in 
the  sediment  as  a  test. 

Construction  of  an  artificially  prepared  ocean  bottom  (kaolin- 
water)  is  recommended  in  order  to  circumvent  the  damage  done  to  the 
cores  when  taken  from  the  bottom  and  to  improve  homogeneity  and 
temperature  stability  of  the  medium. 
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TABLE  I.   Tabulated  Physical  Properties 


SOUND 

SHEAR 

CSED/ 

SAMPLE 

SPEED 
m/sec 

SPEED 
m/sec 

CWAT 

POISSON'S  RATIO 

1. 

1487.0 

22.1 

O.983 

°'^989 

2. 

1502.0 

21.5 

0.993 

°-^989 

3. 

****** 

32.9 

***** 

******* 

4. 

1510.0 

21.1 

0.998 

°-49990 

5. 

1505.5 

23.5 

0.995 

0.4998? 

6. 

1506.1 

26.4 

0.995 

0.49984 

7. 

1517.3 

47.0 

1.000 

o.49951 

8. 

1510.1 

33.5 

0.99S 

o.499?5 

9. 

1493.8 

33.1 

O.987 

o.499?5 

10. 

1520.1 

28.8 

1.000 

0.49982 

******  sound  speed  unattainable 

*  Compr essional  speed  in  sediment/compressional  speed  in  water 
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TABLE    I.      Tabulated  Physical  Properties    (cont.) 


G-l 

G-2 

SAMPLE 

106      a 

dynes/cm 

106    a 

dynes/cm 

DENSITY 
g/cm 

Porosity 

% 

1. 

5.8 

2.5 

1.21 

65.9 

2. 

5.1 

3.0 

1.25 

64.7 

3. 

14.8 

6.9 

1.35 

74.4 

4. 

6.0 

1.8 

1.36 

78.3 

5. 

7.6 

1.1 

1.37 

72.9 

6. 

9.8 

2.7 

1.40 

74.4 

7. 

27.8 

7.9 

I.25 

72.5 

8. 

17.3 

7.0 

1.54 

68.9 

9. 

16.1 

5.2 

1.46 

71.1 

10. 

12.2 

5.8 

1.47 

72.0 
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TABLE    I.      Tabulated  Physical   Properties    (cont.) 


SAMPLE 

PERCENT 

PERCENT 

PERCENT 

NUMBER 

SAND 

SILT 

CLAY 

1. 

2.8% 

29.4% 

67.8% 

2. 

7.2% 

52.6% 

40.4% 

3- 

2.2% 

70.9% 

26.9% 

4. 

6.5% 

25.7% 

67.8% 

5. 

1.4% 

56.5% 

42.1% 

6. 

1.9% 

70.8% 

27.3% 

7. 

2.0% 

72.2% 

25.8% 

8. 

1.9% 

79.6% 

18.5% 

9. 

2.2% 

79.7% 

18 .  l% 

10. 

2.2% 

72.3% 

25.5% 
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TABLE  I.   Tabulated  Physical  Properties  (cont.) 

SAMPLE  pC2 

,,-6  "   2  ,      3     2. 

(10  g-m  /cm  -sec  ) 

1.  2.67 

2  .  2 . 82 

3.  **** 

4.  3.67 

5.  3.10 

6.  3.17 

7.  2.87 

8.  3-52 

9.  3.26 

10.  3.40 

****  Sound  speed  unattainable 
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TABLE    II.      Comparison    of    Shear    Speed  Results 


AVERAGE  STANDARD   DEVIATION      ESTIMATE   OF 

VALUE  OF    DATA   POINTS  PRECISION 


Ratio   of  V./V  0.758  0.062  +    5  % 

is  — 


b 


Shear    wave   speed  29.0   m/sec  8   m/sec  _+   30% 

(from  measured 

rigidities) 

Interface  wave   speed     28.8  m/sec  4  m/sec  _+    10% 

(from    interface   wave 
exper  iment ) 

Shear    wave   speed  33-3  m/sec  4  m/sec  +^  15% 

(from  above) 
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TABLE  III. 
Variation  of  Constants  After  Averaging  at  Low  Viscosities 


^L 

Temperature  I 

:°c) 

4  kHz 

27    kHz 

38    kHz 

10 

6.04 

0.43 

0.34 

17 

5.73 

O.69 

0.62 

22 

6.20 

0.71 

0.70 

K2(K2 

X 

10"3) 

Temperature i 

;°c) 

4kHz 

27    kHz 

38    kHz 

10 

1.36 

1.38 

1.52 

17 

2.96 

2.67 

2.83 

22 

3.81 

2.66 

2.72 
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Figure  1.  Sketch  of  Interface  Wave  Measurement  Apparatus 
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Figure   2.    Visicorder    Prints    from   Interface  Wave 
Measurement 
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Figure  3.  Sketch  of  Viscoelastometer  and  Associated 
Equipment  [2] 
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Figure   4.    Sketch    of   Compress ional   Wave   Velocimeter 
Apparatus    [2] 
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POROSITY  (%) 

Figure   6.    G      as    a    Function    of   Porosity 
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Figure  7-   G,  as  a  Function  of  Sound  Speed 
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Figure   8.      G.    as   a    Function 


of   Poisson's    Ratio 
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Figure   9.    G.,    as    a    Function    of    Density   and  Sound   Speed 
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Figure  12.  G   as  a  Function  of  Density 
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Figure    13.    G„   as    a    Function    of  Porosity 
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Figure    15-    G     as    a    Function    of  Poisson's    Ratio 
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APPENDIX:   CLAY  MINERALOGY 

The  results  of  the  quantitative  analyses  of  the  clay  minerals 
from  X-ray  diffraction  are  shown  in  Table  IV.   The  relatively  high 
percentage  of  chlorite  in  these  samples  is  due  to  river  drainage 
from  areas  of  outcrops  of  the  Franciscan  Formation.   The  graphs  of 
G  and  G      as  a  function  of  the  various  clay  minerals  show  the  data 
from  this  report  as  well  as  those  from  Bieda  (1970) ( Figures  20  -  25) 
Little  in  the  way  of  definite  trends  may  be  seen,  possibly  due  to 
the  fact  that  these  data  are  all  from  nearshore  sediment  samples. 
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TABLE   IV.      Clay   Mineralogy 


PERCENT 

PERCENT 

PERCENT 

PERCENT 

SAMPLE 

MONTMORILLONITE 

ILLITE 

KAOLINITE 

CHLORITE 

2. 

72.88 

13.56 

6.51 

6.98 

3. 

69-31 

13.97 

8.05 

8.66 

4. 

60.72 

17.86 

11.82 

9.60 

5. 

51.35 

21.62 

12.16 

14.86 

6. 

47.33 

25.95 

11.76 

14.96 

8. 

64.91 

12.21 

10.53 

12.37 

9. 

48.82 

24.24 

12.53 

14.41 
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Figure  20.   G   as  a  Function  of  Percent  Montmorillonite 
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